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A nonlocal generalization of Einstein's theory of gravitation is constructed within the framework 
of the translational gauge theory of gravity. In the linear approximation, the nonlocal theory can be 
interpreted as linearized general relativity but in the presence of "dark matter" that can be simply 
expressed as an integral transform of matter. It is shown that this approach can accommodate the 
Tohline-Kuhn treatment of the astrophysical evidence for dark matter. 
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1. Background. In a recent paper the main out- 
lines of a nonlocal theory of special relativity have been 
sketched and some of its experimental consequences have 
been described. The extension of Lorentz invariance to 
accelerated observers in the standard theory of relativ- 
ity is based on a postulate of locality that — on physical 
grounds — is expected to be just an approximation. The 
recognition of the limitations of the locality assumption 
has led to the development of nonlocal special relativity 
(see [l], in and the references cited therein). It is natu- 
ral to consider the extension of this acceleration-induced 
nonlocality to the gravitational field, as the principle of 
equivalence establishes a connection between inertia and 
gravitation. However, a direct approach to a nonlocal 
generalization of general relativity (GR) fails, since Ein- 
stein's principle of equivalence is strictly local. Further- 
more, the heuristic power of Einstein's principle of equiv- 
alence cannot be employed to gain insight into the phys- 
ical connection between nonlocal gravity and nonlocal 
special relativity. This circumstance provides the incen- 
tive to search for other viable approaches to a nonlocal 
generalization of Einstein's theory of gravitation. 

2. GR\\. It turns out that general relativity can also be 
obtained as a special case of the translational gauge the- 
ory of gravity (see [1] ; new developments and further ref- 
erences are contained in recent papers [J as well as mono- 
graphs and ) . This teleparallel equivalent of general 
relativity (GR||) can be succinctly described within the 
framework of a Weitzenbock spacetime. 

In Weitzenbock geometry, the gravitational potential 
is represented by a coframe 1-form -d"' = Ci^dx^, where 
ei°'{x) is the tetrad field. Here ... = 0, 1, 2, 3 denote 
holonomic coordinate indices, while a,/3, ... — 0, 1, 2, 3 
denote anholonomic frame indices. We choose units such 
that c = 1 throughout; moreover, the Minkowski metric 
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tensor rj is given by diag(l, —1, —1, —1). The frame 
e^adi is dual to the coframe; therefore, 
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(1) 



The gravitational field strength is the object of an- 
holonomity 
dx^\ that is, 



holonomity 2-form C" := di?", where C" = iCy"da;* A 
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(2) 



We choose the frame field to be orthonormal throughout. 
The spacetime interval is given by ds^ = gijdx"^ ®dx^ , 
where 



(3) 



A geodesic between two events is defined by the space- 
time path that is an extremum of the spacetime distance 
between them. Finally, in Weitzenbock geometry the 
frame field can be chosen to be globally teleparallel. That 
is, the connection 1-form F^^ = Tia^dx'^ is such that the 
curvature 2-form 
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F„'' A F^'' = ^R,jjdx' A dx^ (4) 



is zero. Therefore, the frames in Weitzenbock spacetime 
can be so chosen that the connection vanishes every- 
where, 



F ^ 



0. 



(5) 



Here the star over the equality sign indicates that the 
equation is valid only in the specially chosen global 
frames. Thus Weitzenbock spacetime has vanishing cur- 
vature but its torsion is in general nonzero and in the 
special global frames reduces to the gravitational field 
strength C^ ". 

3. Gi?|| field equations. The gravitational field equa- 
tions of the translational gauge theory of gravity are 
given by 3 
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where 7^' = \J —g and 

So.' := -\e'o.{C,k^Wp) + Co.k^W'^p ■ (7) 

Here 'hV'-' a are the gravitational excitations that are in 
general linear in the field strengths C'-'q. Moreover, 
the energy-momentum tensor density of matter is given 
by 7^% while has the interpretation of the energy- 
momentum tensor density of the gravitational field. The 
gravitational field equations ^ resemble Maxwell's equa- 
tions. Indeed, the first set of equations expresses the 
divergence of excitation in terms of the sources includ- 
ing gravity (7^* + 8^); that is, all physical fields carry 
energy-momentum. The second set of equations ([5]) is 
"solved" by Eq. ([2]), which expresses the field as the (ex- 
terior) derivative of the potential. 

To recover general relativity in this scheme, we must 
assume that 

= ^ e\ , (8) 

K 

where k = SttG and G is Newton's constant of gravita- 
tion. Here 

^u" 2 ~ + 2e[i"Cj]^^ (9) 

is the modified field strength. It has been shown that 
Eqs. ([H])-® are an equivalent representation of Einstein's 
theory of gravitation Q . 

4. Nonlocal GR\\. The translational gauge theory of 
gravity bears a striking resemblance to electrodynam- 
ics. Just as electromagnetism involves relations con- 
necting excitation = (D, H) with the field strength 
Fij = {E,B), one may look upon Eq. ([5]) as a similar 
constitutive relation in GR||. More generally, one can 
extend constitutive relations to nonlocal ones; in partic- 
ular, it is possible to imagine such relations in vacuum 
in connection with acceleration-induced nonlocal electro- 
dynamics [3] • In a similar way, it is possible to introduce 
nonlocality into the framework of teleparallel gravity. 

Consider, for example, instead of Eq. ([8|), 

xK{x,y)<^,,''{y)V^9¥)d^y], (10) 

where fl{x, y) is Synge's world-function Q, which is half 
the square of the geodesic distance connecting x and y. 
The world-function and related bitensors are all smooth 
functions provided we assume that in the spacctime re- 
gion under consideration each pair of events can be joined 
by a unique geodesic. In Eq. pi?]) . the coordinate indices 
a,6, c, ... refer to x, while i,j,k,... refer to y; moreover, 
we define the bitensors 



and note that ft satisfies the basic partial differential 
equations 

2n ^ g'^'^Qa^b = g'^n.Qj . (12) 

The covariant derivatives at x and y commute for any 
bitensor and it is possible to show that Qaiix,y) — 
^iaix, y) are dimensionless bitensors such that 

lim flaiix, y) = -gai{x) ■ (13) 

Equation (fTO|) introduces the causal scalar kernel 
K[x,y) that indicates the nonlocal deviation from Ein- 
steinian gravity in our model. Here K{x, y) — Q ii x 
lies in the past of y; in Minkowski spacetime, for in- 
stance, this kernel would in general be nonzero only for 
x° > The kernel K is in general a function of the 
coordinate invariants that can be constructed at x and 
y within Wcitzenbock spacctime; for instance, it could 
depend on the Wcitzenbock invariants 

C^jkC''\ Ck.^C'^'' , C,jW'\. (14) 

Other coordinate invariants can be formed, for example, 
from f2 and its covariant derivatives. As an illustration, 
let us consider VL°'ea°'{x) and ri'ei"(2/), which are par- 
ticularly interesting in view of the linear approximation 
scheme that is developed below; indeed, in Minkowski 
spacctime these reduce to x" — y" and y" — x", respec- 
tively. We note that the nonlocal constitutive relation 
pn)) is highly nonlinear in the gravitational variables. Its 
substitution in the field equations ^ then results in a 
nonlocal generalization of Einstein's theory of gravita- 
tion characterized by a given nonlocal kernel K(x,y). It 
should be emphasized that other nonlocal models can be 
considered based on more complicated constitutive rela- 
tions; in fact, Eq. PU)) represents the simplest nonlocal 
constitutive model involving a scalar kernel. The physi- 
cal origin of this kernel is beyond the scope of our con- 
stitutive approach; a more basic theory is needed for its 
determination. 

Let us observe that our approach to nonlocal gravita- 
tion differs essentially from other nonlocal modifications 
of general relativity; see, for example, [l3| and the 
references cited therein. 

5. Linear approximation. The Wcitzenbock spacctime 
in the Ti"'^ = gauge reduces to Minkowski spacetime 
for e.i" ~ Sf. Thus it is possible to develop a linear 
approximation of our nonlocal theory of gravitation by 
assuming 

e," = Sf + ^Z-", , (15) 

where the nonzero components of -0"^ have magnitudes 
that are much smaller than unity. Henceforth ijj^i will 
be taken into account at the first order of approximation, 
where the distinction between holonomic coordinate in- 
dices and anholonomic tetrad indices disappears. The 
metric of the Wcitzenbock spacctime is then given by 

9ij = Vl] + % : = 2V'(ij) , (16) 
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and the gravitational field strengths can now be ex- 
pressed as 



(17) 



+ Sl[^p.^ -^P,\')-Si{t^-i'^\') , (18) 

where = r]ij'ijj'^ . Moreover, the constitutive relation 
(fTOt takes the form 



{x)+nC{x,y)e\{y)cfy, (19) 



where /C(x, y) is the scalar kernel K{x^ y) evaluated in the 
Minkowski spacetime limit. Thus /C is in general nonzero 
for x'^ >y^- The substitution of these relations in Eq. ([6]) 
results in 



dlC{x,y) 
dx^ 



(20) 



This is the main equation of our linearized theory; for 
/C = 0, it reduces to linearized general relativity, since it 
is simple to prove that the Einstein tensor in the general 
linear approximation is given by 



(21) 



We note that diTk = follows in general from the non- 
local field equations The constitutive relation 
is valid for sufficiently weak gravitational fields as the 
constitutive kernel K in the linear approximation be- 
comes completely independent of the gravitational field 
variables. 

6. Reciprocal kernel. Let us now assume that JC{x,y) 
is a function of x — y; then, Eq. (j20p implies that 

G„ (x) + fiCix - 2/)G„ (2/)A = kT,,{x) , (22) 



where we have neglected certain boundary terms, since 
the derivatives of ip'^i are expected to vanish at infinity. 
This is a Fredholm equation of the second kind and it can 
be formally solved by the Liouville-Neumann method of 
successive substitutions. That is, one can move the sec- 
ond term from the left-hand side of Eq. (|22p to the right- 
hand side and then substitute for Gy (y) in the integrand 
its value given by the resulting equation. The repetition 
of this process would lead to an infinite series in terms 
of iterated kernels JCn{x,y), n = 1,2,3,... , that are in 
general defined by ICi{x, y) — IC{x, y) and 



JCn+i{x,y) 



IC{x, z)Kn{z, y)d'^z . 



(23) 



If the resulting infinite series is uniformly convergent [11 
we can define a reciprocal kernel TZ{x, y) given by 



(24) 



such that the formal solution of equation ([22]) can be 
written as 



Gijix) = KTij{x) + K n{x,y)Ti.j{y)d'^y . 



(25) 



To preserve causality, it is convenient to assume that 
the constitutive kernel is of the form 



IC{x - y) - S{x° - 2/0)p(x - y) 



(26) 



It follows from Eqs. (HJ) and (HD) that all of the iter- 
ated kernels as well as the reciprocal kernel are of this 
same form, since the corresponding spatial integration in 
Eq. would be over the entire Euclidean space. In 
particular. 



nix -y)^ 5{x° - y°)qi^ - y) , 



(27) 



where p and q arc reciprocal spatial kernels. 

7. Dark matter. It is possible to express Eq. (j25p as 



G. 



(28) 



and interpret the new source for linearized Einsteinian 
gravity as "dark matter" . Consequently, is the sym- 
metric energy-momentum tensor of "dark matter" that 
is given by 



%jix)= / 7^(x,2/)T,,(2/)A. 



(29) 



Thus "dark matter" is in effect the integral transform 
of matter by the reciprocal kernel TZ(x,y). In our sim- 
ple model, the characteristics of "dark matter" should 
be quite similar to actual matter in accordance with 
Eq. ([^5]) . For example, the "dark matter" associated 
with dust would be pressure-free, while lij is traceless 
for radiation with Tk^ = 0. In particular, for a reciprocal 
kernel of the form (P7|) . wc have for dust of density p, 



PD {t, x)^ J q{x- y)p{t, y)(f 



y, 



(30) 



so that the density of "dark matter" po is in effect the 
convolution of p and q. 

It appears natural to speculate that the astrophysical 
evidence for dark matter may in fact be a manifestation 
of the nonlocal aspect of classical gravity. In the linear 
approximation, the kernel in Eq. I|29p is universal in the 
sense that it is independent of the particular configura- 
tion of matter under consideration. It follows from the 
application of Eqs. and ([M)) to the solar system that 
our predicted dark matter must be a rather small frac- 
tion of matter; otherwise, there would be conflict with the 
solar-system tests of GR. This requirement can be met 
with an appropriate choice of kernel based on the obser- 
vational evidence for dark matter. Consider, for instance, 
the circular motion of stars in the disk of a spiral galaxy 
in connection with the problem of dark matter in such 
galaxies [H, [H, [H, M, M, M, [H, H ■ Outside the 
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bulge, the Newtonian acceleration of gravity for each star 
at radius |x| is toward the galactic center with magnitude 
?;g^/|x|, where vq is the (approximately) constant speed 
of stars in conforrnity with the observed rotation curves 
of spiral galaxies [H, [13, [iBl- The density of dark 
matter in the disk responsible for this behavior is given 
by Vq / (47rG'|xp), which follows from Poisson's equation. 
Neglecting the dimensions of the galactic bulge and set- 
ting p{t,y) ~ MS{y), where M is the effective mass of 
the galaxy, Eq. ([50]) then implies that 



g(x) = 



1 1 

4^|xF 



(31) 



where A = GM/vq is a constant length parameter of 
order 1 kpc. With this reciprocal kernel, the Newtonian 
limit of our linearized nonlocal theory is given by 



AttG 



p(t,x) + 



1 

4^ 



p{t,y)d^y 
|x-y|2 



(32) 



where $ is the Newtonian potential. For a point mass m 
with p{t,x) = m(5(x), Eq. ([5^ has the solution 



Gm Gm ^ /Ixl 



(33) 



where the logarithmic term due to dark matter has a 
negligible influence in the solar system, as required. 

Remarkably, Eqs. (l32l) and ([33|) coincide with the 
Tohline-Kuhn scheme [l9| that, apart from a disagree- 
ment with the empirical Tully-Fisher law, has been quite 
successful in dealing with dark-matter issues in galax- 
ies and clusters — see the lucid review by Bekenstein . 
The universality of kernel (pij) implies M cx Vq-, there- 
fore, the main shortcoming of the linear approximation 
that results in a universal kernel is that the Tully-Fisher 
relation (A/ cx Vq) is violated [20| . 

In a general treatment, an equation of the form ([^5]) 
may still be pertinent, but with a kernel that would 
strongly depend on the particular matter distribution. 
Furthermore, it would be interesting to have a theory 
that would be able to predict the kernel from first prin- 
ciples. But the development of such a theory remains a 
task for the future as it is beyond the scope of the con- 
stitutive approach adopted in this work. 
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